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Abstract 

A derivation of the Boltzmann equation from the Liouville equation by the use of 
the Grad limiting procedure in a finite volume is proposed. We introduce two scales of 
space-time: macro- and microscale and use the BBGKY hierarchy and the functional 
formulation of classical mechanics. According to the functional approach to mechanics, 
a state of a system of particles is formed from the measurements, which are rational 
numbers. Hence, one can speak about the accuracy of the initial probability density 
function in the Liouville equation. We assume that the initial data for the microscopic 
density functions are assigned by the macroscopic one (so, one can say about a kind 
of hierarchy and subordination of the microscale to the macroscale) and derive the 
Boltzmann equation, which leads to the entropy production. 



1 Introduction 

The report concerns the problem of derivation of the kinetic Boltzmann equation from the 
equations of microscopic dynamics (the Liouville equation). The main interest is related 
with the fact that the Liouville equation is reversible in time, while the Boltzmann equation 
does not. On the contrary, for the Boltzmann equation the so-called if-theorem is valid, so, 
it describes the entropy production and time irreversibility. 

The fundamental time irreversibility problem (the Loschmidt's paradox) consists in the 
following: how to reconcile the time-reversible microscopic dynamics and the irreversible 
macroscopic one. One of the problems in the framework of the irreversibility problem is the 
derivation of the Boltzmann equation from the Liouville equation. 
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An elegant derivation of the Boltzmann equation from the Liouville equation has been 
proposed by Bogolyubov []]. He uses the BBGKY (Bogolyubov-Born-Green-Kirkwood- 
Yvon) hierarchy of equations, thermodynamic limit and some additional assumptions (which 
do not follow from the Liouville equation). This derivation leads to the divergences in high 
order corrections to the Boltzmann equation [2J. 

Another derivation of the Boltzmann equation was proposed by Lanford [5] . He also uses 
the BBGKY hierarchy of equations, the Grad (or Boltzmann-Grad) limit, all assumptions 
are included in the initial conditions. But the derivation can be applied only to small times. 
On the contrary, the Boltzmann equation is interesting from the viewpoint of the large time 
asymptotics (the relaxation to the Maxwell distribution). 

We propose a new derivation of the Boltzmann equation. The BBGKY hierarchy of 
equations are supplemented with the ideas of the functional mechanics, recently proposed 
by I. V. Volovich [HE]. Also we introduce two scales of space-time: a macro- and a microscale 
(a kinetic scale). 

A kinetic equation for the system of two particles has been obtained in [B]. But the 
obtained equation is time-reversible, as the Liouville equation. In the present work we 
obtain the irreversible Boltzmann equation, which leads to the entropy production. 



2 Liouville and Boltzmann equations 

Let N particles in a region GcR 3 with the volume V be given. Their state in an arbitrary 
moment of time t is described by a function f(xi,x 2 , ■ ■ ■ ,x^,t), where %i = (qi,Pi), 
(a position), pi G R 3 (a momentum), i = 1,2, ... ,N. Let the normalization condition be 
J f dxi . . . dxjy = V N , i.e., y^f is a probability density function of the iV-particle system. 
The dynamics of the function / is given by the Liouville equation 

f = W/}, (i) 

where 

N 2 N I I N 

H ^p) = E£ + E + E ( 2 ) 

i>j 

is a Hamiltonian of the system, m > is the mass of a single particle, {-, -} are the Poisson 
brackets. The first term in H corresponds to the kinetic energy of the system, $(r) is the 
interaction potential of the particles, U{qi) is an external potential. 

We assume that the function $(r) is continuously different iable, bounded from below, 
$(r) — > +oo asr — > 0, $(r) — > asr — > oo. Also for simplicity we will assume that 
$(r) monotonically decreases as r increases (which corresponds to a repelling force). Let the 
function U be continuously differentiable, U (q) — > +oo as q — > dG, where dG is the boundary 
of the region G (so, the external potential does not allow the particles to go beyond G). 

Further, \x > is a small dimensionless parameter. It has a sense of the relation of 
the interaction radius ro to the mean free path Z, i.e., /i = j < 1. This parameter in the 
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Hamiltonian means that the particles interact on the scale much smaller than the scale of 
kinetic events. 

Time reversibility of the Liouville equation is expressed in terms of the following 

Proposition 1. Let f(q,p,t) be a solution of the Liouville equation (Tjp. Let the Hamiltonian 
satisfy the equality H(q,p) = H(q,—p) (for example, the Hamiltonian has the form (2)). 
Then f(q, —p, —t) is also a solution of the Liouville equation. 

Note, however, that in spite of the formal reversibility, the solutions of the Liouville 
equation obey the property of delocalization which can be regarded as irreversibility [U [5] . 

Consider now a gas from the point of view of kinetic theory. In this case the gas is 
described in terms of a single-particle density function f(q,p, t). Let it be normalized on the 
volume V (i.e., f(q,p,t)/V is a probability density function), q 6 G, p e R 3 . One of the 
fundamental equations of kinetic theory is the Boltzmann equation: 

?l = _L?l + audf +Si . f (3) 

dt m dq dq dp 
f \P — Pi\ 

St/ = n / [f(q,P,t)f(q,p' 1 ,t)-f(q,p,t)f(q,p 1 ,t)]dadp 1 , 

jr 2 xr3 rn 

where n > is the mean concentration of the particles (the mean number of particles in a 
unit of volume, da = rdrdip. Also p' and p[ are the momenta that will have two particles long 
after the collision provided that they had the momenta p and p\ long before the collision 
with the impact parameter of the collision r and the polar angle (p. So, (r, (p) e IR 2 are 
polar coordinates on the plane perpendicular to the relative velocity vector (p — pi)/m. 
Thus, p' = p'(p,pi,r,ip), p[ = p[(p,pi,r,ip), the dependence is defined by the two-particle 
Hamiltonian without the external potential 



2m 2m ' fi 



flg = il. + ftL + g ( |gi-ftl ). ( 4 ) 



The expression St / is called the collision integral. 

This is an important nonlinear equation which describes the relaxation of the function 
/ to the Maxwell distribution [Tj. One of the properties which can be easily proved is the 
so-called Boltzmann i?-theorem: 

Proposition 2 (if-theorem). Let f(q,p,t) be a solution of the Boltzmann equation |5|) and 
the quantity 

f(q,P,t)kif(q,p,t)dqdp, 



r 



Qy = G x M 3 ; be well-defined (i.e., the integral converges). Then > 0. 

This proposition states the entropy production (the quantity S(t) = — H(t) can be re- 
garded as the entropy of the gas) and, hence, the irreversible character of the gas dynamics. 

Thus, we have obtained two contradictory conclusions: if we consider the g whole 
(in terms of a single-particle distribution and the Boltzmann equation), the dynamics is 
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irreversible in time, while if we consider the gas as a system of a finite number N of particles 
(in terms of an iV-particle distribution function and the Liouville equation), the dynamics 
is reversible. 

It seems that the reductionism does not work, the time irreversibility is a property of 
macrosystems which cannot be reduced to the microscopic level. So, the problem about 
another type of relation between the different levels of description arises. 

Let us return to the Liouville equation. According to Bogolyubov, let us define the 
s-particle distribution functions f s , s — 1,2, . . . , N — 1: 



where fiy = G xl 3 is the phase space of a single particle, /jv = /• Usually, the function / in 
the Boltzmann equation is associated with the single-particle function fx. This is a common 
point in the derivations of the Boltzmann equation according to Bogolyubov, Lanford and 
others. We also follow this idea, but supplement it with the idea of subordination of different 
space-time scales. 

3 Micro- and macroscale 

Both Bogolyubov and Lanford start with the Cauchy problem for the Liouville equation: 



But we should keep in mind that the initial distribution function f° is not given "objectively" , 
but constructed based on measurement results. Note that the measurement results are 
rational numbers (this is a starting point of p-adic mathematical physics [EJ|9]). See [10] for 
the detailed description of the construction of the probability density function starting from 
the rational-valued measurement results. See also [11] for the discussion of the functional 
dynamics of a system under often measurements. 

The accuracy of the measurements is essential here. If our measurement instruments 
allow us to register variations of the physical values (and, hence, the initial distribution 
function /°) on the scale of particle interaction radius, the solution of the Cauchy problem 
has a physical meaning. 

But in practice, if we consider the kinetic events, our measurement instruments can 
register variations of the initial distribution function only on the scales much larger than 
the interaction radius [H [12] ■ The scale of order of the interaction radius r will be called 
"microscale" , the scale of order of the mean free path I - "macroscale" (or "kinetic scale" ) . 
fi = y <C 1 is a scaling parameter. Our measurement instruments can register the variations 
of the physical values only on the kinetic scale. 

In this case we do not know the iV-particle distribution function f°(x\, . . . ,xn), since 
it reflects the information about correlations of the particles on the scale of order of the 
interaction radius, which cannot be registered. This information is essential. For example, 





(5) 
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if the distance between the particles is microscopically large (|<ft — Qj\ 3> tq for all i 7^ j), 
there is no correlation between them and the following factorization property is satisfied: 

N 

f°(x 1 ,x 2 ,...,x s ,0)=]Jf°(x i ,0). 

Exactly this property will be registered by the instrument, since it can register only the 
macroscopic variations. However, if we extrapolate this property over the whole phase space, 
we obtain an infinite mean energy E = f u Hf dx\ . . .dxisf (since $(r) — > +00 as r — > 0). 

So, the correlations between the particles are essential on the microscale, but the in- 
strument cannot register them. Hence, we do not know the initial iV-particle distribution 
function f° and the Cauchy problem §5§ has not a direct physical meaning. 



4 The main theorem 

Let us formulate another problem for the Liouville equation which will be in accordance 
with the physical meaning. The instrument allows us to establish the initial single-particle 
distribution function f®(xi). One can say that /° is a "macroscopic" probability distribution, 
since it varies on the scale much larger than the particle interaction radius. Let /° be 
continuous with its partial derivatives over each argument. 

Let us define the following problem for the Liouville equation: 



S-ittWxi, x 2 , t - At) - h„{x u t - At)f ltl (x 2 , t - At)} ->• 

(as ix — > 0, iV — > 00, iV/i 2 = const, and At — > 0, — 

fi 

hn{x x , 0) -> fl(x x ) (as n -> 0), 



At — > 00) 



(6) 



O-V 



Hf^ dx\ . . . dxN < 00, 



k f li (xi,...,x N ,t) = f^x l{1) ,...,Xi iN) ,t), (i(l),...,z(JV)) = P(1,...,N), 

Here, is the two-particle Hamiltonian flow, i.e., S^ 2 ' ip(xi,x 2 ) = l p(xit,x 2 t), where 
{xu, x 2t ) is the phase point (xi, x 2 ) moved along the flow defined by the Hamiltonian © 
on t, -P(l, . . • , N) is a permutation of the numbers 1, . . . , N. The second condition in ([6]) is 
understood as a weak limit [13] over the variable p 2 , i.e., 



lim / S {2) At [f 2fl {x 1 ,x 2 , t- At) - h^x^t - At)f lfl (x 2 , t - At)}<p{p 2 )dp 2 = 

for any function cp such that the integral exists, f ^ depends on /x, since /x is a parameter in 
the Hamiltonian. 

The limit /j — >■ 0, iV — > 00, iV/i 2 = const is called the Grad or Boltzmann-Grad limit [Hj 
and was used by Lanford. 
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Theorem. Let the function f^(x\,. . . ,x^,t) satisfy the problem and fi^(xi,t) tend to 
some function fi(xi,t) in the Grad limit in every point x± G Qy Let $(r) be monotonically 
decreasing function and lim r 7 $(r) = C 7^ 0, 7 > 2. Then the function fi(x\,t) satisfies 

r— >oo 

the Boltzmann equation: 

W± = _L^ll + ^L^ll + st fl 

dt m dq dq dp 
Stf = nfi 2 / [fi{q,p',t)fi{q,p[,t) - f 1 (q,p,t)f 1 (q,p 1 ,t))dcrdpi, 

JK 2 xR 3 m 

f(x 1 ,0) = f°(x 1 ). Heren = ^. 

Hence, the if -theorem and the entropy production (see Proposition 2) are also valid. 

Let us discuss the formulation of problem ([6]). The third condition is the initial data for 
the single-particle function f\. The fourth condition means the fmiteness of the energy, the 
fifth condition means the symmetry of the density function with respect to permutations 
(i.e., the particles are indistinguishable). 

The most interesting and crucial is the second condition. The knowledge of the initial 
single-particle function fi(x\) is not sufficient to get a unique solution for the single-particle 
function fi(xi, t). Since the condition of the form f(xi, . . . , xn, 0) = f°(xi, . . . , x^) has not 
a direct physical meaning, we must have some additional condition. 

As we said above, there are two scales of space-time in this consideration: the microscale, 
related to the interaction of the particles, and the macroscale, related to the kinetic phenom- 
ena. The two-particle function f 2 relates to the microscopic scale, because it incorporates the 
information about the pairwise correlations of the particles on the distances of order of tq. 
The single-particle function f\ relates to the macroscopic scale, since it does not incorporate 
the information about the particles' correlations and the kinetic theory is expressed in terms 
of this function. The second condition in fl6]) means that the initial value for the microscopic 
function f% are assigned by the macroscopic function f\. So, instead of the specification of 
the initial microscopic function /°, we specify only the initial macroscopic function f® and 
impose a condition on the microscopic function: in a certain sense it is subordinated to the 
macroscopic one (in sense that its initial values are assigned from the macroscale). 

Let (q,t) be macroscopic space-time variables, Aq ~ I, At ~ l/u, where u is the mean 
velocity of the particles, and let (£, r) be microscopic space-time variables, A£ ~ r , At ~ 
Tq/u. These variable are related to each other by the scale transformation 

p, /i I 

With such a transformation a macroscopically infinitesimal region may be infinitely large 
from the microscopic point of view. Exactly this situation we can see in the second condition 
in (JHJ): At — > 0, but At — — — > 00. Of course, these two scales meet each other in 
the collision integral of the Boltzmann equation: a collision is considered as a point and 
momentary act on the macroscale, but it takes place on the infinite space during the infinite 
time on the microscale. 

The great disparity between microscopic and macroscopic scales as one of the origins of 
the irreversible macroscopic behaviour was pointed out in (T5l [16]. The used rescaling of 
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the space-time ([7]) is typical for the derivation of the Boltzmann equation from the Liouville 
equation [T7J [18] . In the case of lattice dynamics this rescaling was used in [191 120] for the 
derivation of kinetic and hydrodynamic-type equations. Our proposition is to introduce the 
subordination of different space-time structures expressed in the form of the second condition 
in OH]), which gives a new way of derivation of the Boltzmann equation. 

One can say about a kind of hierarchy and subordination of the scales: the initial values 
for the processes on the microscale (interactions of the particles) are assigned from the 
processes on the macroscale (kinetic phenomena). Note that the idea of the hierarchy of 
times (namely, the microscopic, kinetic and hydrodynamic relaxation times) in a slightly 
different sense was first proposed by Bogolyubov pQ. 

Remark. We can note that in our limiting process the overall mass of the gas M = mN tends 
to infinity, since m is constant and N tends to infinity. However, we can also rescale the 
mass of a single particle in the way m M = fi 2 m and rescale the interaction potential in the 
way 3>/i(r) = yU 2 $(r) as it was proposed by Grad [14J (we substitute m by m M and $(l^^i) 

by $ /Lt ( ) in ()2])). The Hamiltonian system with these additional rescalings can be easily 
reduced to the considered one. So, the gas is considered as an infinite number of negligible 
mass particles in the limit. This exactly corresponds to the macroscopic intuition. 

Finally, we would like to note that the functional approach can be useful in the method of 
molecular dynamics simulations. At present the molecular dynamics follows the Newtonian 
approach and simulates the movement of material points. On this way it is hard to obtain 
the properties of complex and, moreover, biological systems, which constitute the aim of 
the molecular dynamics. Furthermore, the problem of uncontrolled cumulative errors in 
numerical integration is known. The simulation of the motion in terms of the probability 
density function seems to be more appropriate for obtaining the properties of the complex 
and biological systems and more justified from the computational point of view. Also the 
present report suggests that the initial conditions for the probability density function should 
be chosen not in an arbitrary way, but rather in the form suggested by the system as a whole. 
In our case such conditions are given by (EJ). Let us note that these conditions are statistical 
in essence, they are formulated for the probability distribution and cannot be reduced to the 
conditions on the initial positions and momenta of the individual particles. 

5 Proof of the theorem 

The proof of the theorem is divided into four lemmas. Expressions like A pa B means 
A — B — > (in the Grad limit iV — > oo, fi — y 0, Nfi 2 = const). Also for simplicity we will 
skip the lower index /i of the functions etc. So, everywhere in this section f s means 

fs/it s 1,2,.... 

In fact, the presented theorem and its proof is a rigorous mathematical formulation of 
the variant of the Bogolyubov's derivation presented in [T2] . The distinguishing of two scales 
of space-time and the Grad limit allow to do this. 
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Lemma 1. The functions f s (xi,t), s = 1,2, . . . ,N — 1, satisfy the equations 

— = {h s , fs] + — — / { 2^ n — - — ), fs+i} dx s+l . 

J S7t/ a — i i 



L v i=l 

Here 

i=l i,j=l i=l 

is t/ie s-particle Hamiltonian. 

In particular, the function fi satisfies the equation 



dfr p 1 df 1 , dU{ qi )dh , N-l f d^C-3^) d f 2 



dt m dqi dq\ dpi V J n dq\ dpi 



(xi, x 2 , t) dx 2 . 



The proof is straightforward: integration of the Liouville equation over the variables 
x s+ i, . . . , xn- This well-known set of equations is called the BBGKY hierarchy. 

Lemma 2. In the sense of weak limit the following limiting equality is satisfied: 

f 2 (xi,x 2 ,t)nf l (X 1 ,t)f 1 (X 2 ,t), 

where X { = {Q h Pi) = (q i0 + ^ , p i0 ) , % = 1,2, (x 10 ,x 20 ) = S { X( Xl ,x 2 ), At^0,f^oo. 
Here x i0 = {q io ,Pio), z = 1,2. 

Proof. Firstly, we want to proof the equality 

f 2 (x 1 ,x 2 ,t) « fi(x w ,t- At)/i(a;2o,*- At), (9) 
Indeed, according to lemma [fl the function f 2 satisfies the equation 



3- 



8=1 " U V j=i 

Hence, 

f 2 {xi,x 2 ,t) = S m At f 2 (x 1 ,x 2 ,t - At)+ 

+ f fEW./^^T)} + ^/ {^$(^ :i ^),/3(xi,x 2 ,a; 3 ,r)}dx3 

Since At — >■ 0, we can neglect the integral term: 

f 2 {xi,x 2 ,t) « S^J^x^x^t - At) = f 2 (x l0 ,x 20 ,t - At). 
By the second condition in (JB]), 

^(ario^ao,*- At) « /i(xi ,t- At)f\(x 20 ,t- At). 

8 



dr. 



Equality © has been proved. 

Now we apply these arguments once again: 

h(qio,Pm,t- At) = ftfan + ^—^,p i0 ,t)- 

m 



{U(q 1 ),fi(x i0 ,T)}+ N 2 [ {$( — —),f2(XiQ,X 3 ,T)}dx 3 



t—At 



V Jn v V 



dr. 



□ 



Since At — > 0, we can neglect the integral term: 

fi(Qio,Pio,t - At) « /a(g i0 + ^—^,Pio,t), 

m 

so, 

f 2 (x 1 ,x 2 ,t) & f 2 (x 10 ,x 20 ,t- At) w/i(Xi,t)/i(X 2) t). 



Corollary. 

, 5 $(fa^) a/2 , 3$(]»=2ll) 9 

/ k-^ —{x ll x 2l t)dx 2 K> I — — [/i(Xi,t)/i(X 2 ,t)]dx 2 , 

JfV c?i a Pi Jn v OQi a Pi 

and 

dfi„ Pidh du{ qi ) dh f d<$>{^) o 

-kt~ 7^— + ^ TT + n a ■^-[fi( x u t )fi( x 2,t)\dx 2 10 

dt mdqi dqi dpi J Uv dqi dpi 

(n = y , we have replaced the factor N — I by N, because the integral has an order of jj 2 and 
tends to zero). 

Equation ({TO]) is the so-called Bogolyubov equation. Bogolyubov claims that this equa- 
tion is more precise than the Boltzmann one and starts with this equation (not from the 
Boltzmann equation) in order to derive the hydrodynamic equations with viscosity and heat 
conduction. 

Lemma 3. 

~' ' >><> ' l> (>l ( ' /i>( ^ 1 +n I — — -S-\fi{qi,Pio,t)fi(qi,p3o,t)]dx2, 



dt m dqi dqi dpi J n dqi dpi 
where fl = K 3 xR 3 . 

Proof. Since At — > 0, we have Qi m q^ % = 1,2. The function fi tends to a continuous 
function as p, -» 0, so, fi(Qi,p i0 ,t) « fi(qi,Pio,t). 
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Then, replace the integration domain fly by Jl r ((?i) = B r (qi) x R 3 3 (q 2 ,p 2 ), where 
B r (qi) C R 3 is a ball with the center in qi and the radius r, r — )■ 0, r — > 00. We can do such 
replacement because of the asymptotic properties of the potential $: 



/ 

Jn 



[fi{qi,Pw,t)fi{q 2 ,p 20 ,t)}dx 2 



dqi dpi 



J 

Jn 



-[fi(qi,Pw,t)fi(q 2 ,P20,t)}dx 2 . 



Slr(qi) 9 Ql d Pl 



Now q 2 — q\ — > in the integration domain, since r — > 0. So, we can replace fi(q 2 ,p 2 o,t) by 
fi(qi,p 2 o,t). Finally, due to the properties of $ we can once again replace the integration 
domain Q r (qi) by Q. □ 

So, we have replaced the spatial argument of the functions fx by the constant q\. Note 
that the integrand still depends on q 2 , because p i0 = Pio(qi,Pi,q 2 ,P 2 )- Note also that the 
mathematical condition that f\ is a continuous function in the limit // — > means from 
the physical point of view that the single-particle distribution function varies significantly 
only on the kinetic scale (the variations on the scale of order of the interaction radius are 
negligible) . 

Lemma 4. 

n / ^ ^— [fi(qi,Pio,t)fi(qi,P 2 o,t)} dx 2 « St/i 

</fV opi 

Proof. Since gi is a constant in the terms f\{qi,Pio, t), i = 1,2, let us introduce the function 
g(p,t) = fi(qi,p,t) and rewrite the integral in the form 



/ ^ -k— [fi(qi,Pw,t)fi(qi,p 20 ,t)]dx 2 = / — [y(pio,05(P20,0]^2- 

./n dqi dp! J n dqi dp! 



The two-particle Hamiltonian without the external potential can be represented in two 
ways: 

2m 2m // 2m 2m \x 2m 2m 

The last limiting equality is satisfied, since ^ — > 00 and, hence, ^ 20 ~ 910 — > 00. Physically 
this is just the law of conservation of energy: at the moment t two particles obey the kinetic 
energy and the potential energy of interaction, at the moment t — At they are far from each 
other (relative to their interaction radius), so, they have only the kinetic energy. 
Now, from the one side, 

2 2 

{Hi g(p 10 ,t)g(p 20 , t)} » + ^,g(p w ,t)g(p 20 ,t)} = 0. 
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From the other side, pio and p 2 o are functions of qi,pi, q 2 ,p 2 , and 



0$(l22z2li) q d 

{H 2 , g(pio, t)g(j) 20 , t)}= 7^ — ^— [g{Pw, t)g(p 2 o, t)]+ jr^ — ^—\g(pw, t)g(p2o, t)]- 

oqi opi oq 2 op 2 

-K-[g(pw, t)g(p2o, *)] ^— [gipw, t)g(pw, *)]• 

m oqi m dq 2 

So, 

7^ — ^—[g(pw,t)g(p 20 ,t)) + 7—^ — — \g(pio,t)g(p 20 ,t)]- 

oqi dpi dq 2 dp 2 

- —^-[g(Pio,t)g(p 20 ,t)] - —-^-[g(pio,t)g(p 20 ,t)] ~ 0. 
m oqi m dq 2 

Express the first term from the others and substitute it to the integral: 

d r , jw MJ f P2 - Pl d 



dqi dp i 



[g(pw,t)g(p 20 ,t)]dx 2 = j P2 — Pl [g(pw, t)g(p 20 , t)]dqdp 2} 



where q = q 2 — q\ (the term with the derivative ^ vanishes after the transformation of 
the corresponding integral into the integral over the surface in the momentum space, we use 
the integrability of the function g(p,t) in p). Let us introduce the cylindric coordinates q = 
(z, r, if), where the axis z is directed along the vector p 2 —pi- Then (p 2 ~Pi)-§^ = \p 2 — Pi\-§^- 
By integration over z we get 

n f— — —Tr[g(Pio,t)g(p 20 ,t)}dqdp 2 = n [ — — — [g(pi , t)g(p 20 , t)] 
Jn m dq J D m 



rdrdipdp 2 

-oo 

n/i 2 / [g(p[, t)g(p 2 , t) - g(p h t)g(p 2 , t)}dadp 2 = St f u 

J D m 



where 

Pi = lim pi (xi,x 2 ), p[= lim p i0 (xi,x 2 ), 

z— >— 00 z— >+oo 

i = 1,2, are the momenta before and after the collision correspondingly, da = pdpdip is a 
differential cross-section of the collision, p = -, D = M 2 x R 3 . 

This concludes the proofs of the lemma and the theorem. □ 



6 The problem of the existence of solution 

The main problem of the presented derivation is the following: whether exists a solution of 
problem ([6]) such that fip(xi,t) tends to a continuous function fi(xi,t) in the Grad limit 
for every moment of time t (this is a condition of the theorem). 

This condition is analogous to the Boltzmann's molecular chaos (or molecular disorder) 
hypothesis and to the Bogolyubov's condition of the absence of correlations between two 
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particles before their collision. In fact, we ask whether exists a solution of the Liouville 
equation such that the most crucial second condition in ([6]) is satisfied at every moment of 
time. Or, equivalently, whether exists such family of initial functions f®(x\, . . . ,Xn), \i > 0, 
such that fi tends to the predefined function f® as /x — > and the second condition in <^ 
is satisfied at every moment of time t > 0. Intuitively it seems true, but the problem is to 
prove it rigorously. 

A rigorous proof of this statement could be done by the representation of a solution of 
the BBGKY hierarchy in the form of the Dyson expansion series. This was used by Lanford 
[3]. The Dyson expansion series supports some intuition: every term of the expansion is 
related to the number of collisions. Also in some sense the Dyson series is related to the 
path integral. But unfortunately, the Dyson series converges only on small times. For this 
reason the Lanford's derivation of the Boltzmann equation also applies only to small times 
(see the Introduction). 

The fact that the Dyson series can be applied only to small times (less than the mean free 
time) because of secular terms was noticed also by Bogolyubov pQ. In order to overcome this 
problem he proposed a more sophisticated method based on some additional assumptions 
and methods of nonlinear dynamics. But there is no idea how these additional assumptions 
could be proved. Besides this, as we said in the Introduction, his method leads to the 
divergences in the high order corrections to the Boltzmann equation. 

Note that the limit in the second condition in (jH]) is formulated in the sense of weak limit. 
This limit is not satisfied in the pointwise sense for an arbitrary moment of time. But there 
are arguments that this limit can be satisfied in the weak sense, which is sufficient for us 

01221. 

Note also that we assume the pointwise convergence of the function /^(xi, t) to the 
function f\{x\,€) as \x — > 0. However, the theorem and the proof are still valid if this 
convergence also is valid only in the week sense. 

7 Conclusions 

In this report we have proposed a derivation of the Boltzmann equation from the Liouville 
equation based on some ideas of functional mechanics and measurement theory. The main 
features of this derivation is the distinguishing of two scales of space-time (micro- and macro- 
scopic or kinetic) and the subordination of the processes on the microscale to the processes 
on the macroscale. 

According to the traditional paradigm of the mathematical physics, the dynamics is 
completely determined if we know the law of motion, i.e., a differential equation, and the 
initial values for it. However, the initial values themselves are understood as something 
external to the equations of mathematical physics ("As regards the present state of the 
world. . .the laws of nature are entirely silent" [21] ). We propose another picture: the initial 
values for a given level of nature are assigned from the higher level. 

So, instead of the reductionism, which claims the reducibility of all levels of the nature to 
the most microscopic level, we propose a hierarchical picture of the world: the lower levels 
of the nature are subordinated to the higher ones. 
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